STM/STS Study on 4a x 4a Electronic Charge Order and Inhomogeneous Pairing Gap 

in Superconducting E^S^CaC^Og+a 
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We performed STM/STS measurements on underdoped Bi2212 crystals with doping levels p ~ 
0.11, ~ 0.13 and ~ 0.14 to examine the nature of the nondispersive 4a x 4a charge order in the 
superconducting state at T <C T c . The charge order appears conspicuously within the pairing gap, 
and low doping tends to favor the charge order. We point out the possibility that the 4a x 4a charge 
order will be dynamical in itself, and pinned down over regions with effective pinning centers. The 
pinned 4a x 4a charge order is closely related to the spatially inhomogeneous pairing gap structure, 
which has often been reported in STS measurements on high-T c cuprates. 

PACS numbers: 74.25.Jb, 74.50.+r, 74.72.Hs 



INTRODUCTION 



Clarification of the nature of the pseudogap state is ex- 
pected to provide a clue to understanding the mechanism 
of high-r c superconductivity. The pseudogap state ap- 
pears even in lightly doped regions of Ca2- x Na x Cu02Cl2 
(Na-CCOC) and Bi 2 Sr2CaCu20 8+<5 (Bi2212), where the 
pseudogap is of an asymmetric, V-shaped type at very 
low temperatures and has been referred to as the zero 
temperature pseudogap (ZTPG). 1,2 It was recently re- 
vealed by STM/STS studies that in the ZTPG regime a 
nondispersive 4a x 4a charge order appears in the two- 
dimensional (2-d) spatial map of energy-resolved differen- 
tial tunneling conductance dl/dV, which is proportional 
to the local density of states (LDOS). 1 ' 2 A nondispersive 
~ 4a x 4a charge order, electronic in origin, was also re- 
ported in the LDOS maps measured in the pseudogap 
state of Bi2212 at T > T c . 3 Such a spatial structure in 
the LDOS maps was first observed around the vortex 
cores of Bi2212 exhibiting pseudogap-like V-shaped STS 
spectra with no peaks at the gap edge. 4 ' 5 Such charge 
orders have attracted much attention because the charge 
order can be a possible electronic hidden order in the 
pseudogap state. 1,3 

In measurements of LDOS maps in the superconduct- 
ing (SC) state of Bi2212, Hoffman et al. and McElroy 
et al. found a strongly dispersive 2-d spatial struc- 
ture, which has been successfully explained in terms of 
SC quasiparticle scattering interference. 6-8 Furthermore, 
Howald et al. and Fang et al. reported a nondisper- 
sive ~ 4a x 4a charge order with anisotropy in the SC 
state of Bi2212 in addition to weakly dispersive ones, and 
claimed that the charge order was due to the stripe order 
and coexisted with the superconductivity. 9,10 However, 
the nondispersive ~ 4a x 4a charge order at T < T c was 
not confirmed in other group's LDOS measurements on 
Bi2212. 3,6 Very recently the nondispersive 4a x 4a charge 
order at T < T c was found to appear in heavily under- 
doped SC Bi2212. n The charge order is commensurate 
(4a x 4a) and has an internal structure with a period 
of 4a/3 x 4a/3, which is just like the electronic charge 



order reported by Hanaguri et al. for lightly doped Na- 
CCOC. 1 The observation of almost the same charge order 
for both cuprates Na-CCOC and Bi2212 provides definite 
evidence that the nondispersive 4a x 4a charge order de- 
velops on the Cu-0 layer. The 4a x 4a charge order is 
likely to be dynamical in itself, and pinned down over 
regions with effective pinning centers. 11 To understand 
the nature of the nondispersive 4a x 4a charge order, it is 
desirable to investigate the charge order on crystals with 
different doping levels and/or pinning centers of different 
natures. 

In the present work, we studied the nondispersive 
4a x 4a charge order from STM measurements on Bi2212 
crystals with different doping levels and/or pinning cen- 
ters of different properties. We found that low doping 
tends to favor the development of the 4a x 4a charge order 
though it would be dynamical without pinning centers. 
We also studied the STS spectra over the same region 
where STM images were taken, and point out that the 
spatially inhomogeneous gap structure, often reported in 
STS measurements on high-T c cuprates, 12-18 will corre- 
late with the appearance of the pinned ~ 4a x 4a charge 
order. 



II. EXPERIMENTAL PROCEDURES 

In the present study, single crystals of Bi2212 were 
grown using traveling solvent floating zone method. We 
estimated doping level p of the Cu-0 layer from the SC 
critical temperature T c determined from the SC diamag- 
netism and the characteristic temperature T max of the 
normal-state magnetic susceptibility; both T c and T max 
follow empirical functions of p. 19,20 The doping level was 
controlled by changing the pressure of oxygen atmosphere 
in the course of growing the crystal. We performed 
STM/STS measurements on three different single crys- 
tals a (p ~ 0.11, T c ~ 72 K), (3 (p ~ 0.13, T c ~ 78 K) 
and 7 (p ~ 0.14, T c ~ 81 K), and report the results on 
typical sample pairs (A, B), (C, D) and (E, F) cut from 
a, P and 7 single crystals, respectively. In the present 
STM/STS experiments, Bi2212 crystals were cleaved in 
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an ultrahigh vacuum at ~ 9 K just before the approach 
of the STM tip toward the cleaved surface in situ. Bi2212 
crystals are usually cleaved between the upper and lower 
layers of the Bi-0 bilayer. In Bi2212 crystals, excess oxy- 
gen atoms contained within Bi-0 bilayers provide Cu-0 
layers nearby with mobile holes. However, excess oxy- 
gen atoms will be appreciably lost during the process of 
cleaving the crystal at high temperatures. So, to sup- 
press the loss of excess oxygen atoms, i.e. mobile holes, 
to as low a level as possible, we cleaved the crystals at low 
temperatures (~ 9 K). In the present study, STM images 
(512 x 512 pixels) were measured over the surface areas 
of ~ 38 nm x 38 nm for samples A, C, D, E, ~ 23 nm x 
23 nm for sample B and ~ 19 nm x 19 nm for sample F 
in the constant height mode under constant sample-bias 
voltage V s applied between the tip and the sample. We 
were able to observe atomically resolved STM images at 
various bias voltages from a low bias of V s = 10 mV to a 
high bias of 800 mV. The differential conductance dl/dV 
was measured by using a standard lock-in technique with 
an ac bias modulation of 3 mV and a frequency of 4 kHz. 



III. RESULTS AND DISCUSSION 

A. STM images of Cu-O layer; the 4a X 4a charge 
order 

The cleaved Bi-0 layer of Bi2212 crystals is semicon- 
ducting, with a gap of the order of 0.1 eV (Abi-o)- 
Therefore, if we choose a high bias voltage V s , which lies 
outside the semiconducting gap Aei_o where the elec- 
tronic states exist in Bi-0 layers, in the STM experi- 
ment, the STM electron-tunneling occurs predominantly 
between the STM tip and the cleaved Bi-0 layer (Fig. 
1(a)). Thus we can observe the Bi-0 layer selectively in 
STM imaging at a high bias (V s > Aei-o/e), when we 
keep the STM tip at a distance from the sample surface. 
On the other hand, if we choose a low bias V s , which 
lies within the semiconducting gap Aei-o where elec- 
tronic states do not exist in Bi-0 layers but in the Cu-0 
layer , the STM electron-tunneling can occur between the 
STM tip and the Cu-0 layer which is buried just below 
the cleaved Bi-0 layer (Fig. 1(b)). In STM imaging at 
a low bias (V s < Asi-o/e), we can observe the Cu-0 
layer selectively when the STM tip approaches the sam- 
ple surface so that wave functions of carriers between the 
STM tip and the Cu-0 layer can overlap. 21 In fact, STM 
images taken on sample A at high and low biases had 
different features, especially, with respect to the missing 
atom rows inherent in the Bi-0 layer, as shown in Fig. 
2(a). 22 In the STM image taken at 600 mV (high bias), 
the missing atom rows appear very clearly (the inset of 
Fig. 2(a)). On the other hand, the missing atom rows be- 
come very weak in the STM image taken at 30 mV (low 
bias), as seen in Fig. 2(a). The latter result confirms that 
in the low-bias STM imaging the STM tunneling mainly 
occurs between the STM tip and the Cu-0 layer with no 



missing atom rows. 




low-bias 

FIG. 1: (Color online) Schematic energy diagram of Bi2212 
and illustration of STM measurements at (a) high (V s > 
Abi-o) and (b) low biases (V s < Aei-o). The density of 
states N(E) for Cu-O and Bi-O layers are schematically rep- 
resented by thin and thick lines, respectively. In the high-bias 
STM experiment, electron tunneling occurs predominantly 
between the STM tip and the Bi-0 layer when the tip-sample 
separation is large, whereas in the low-bias STM experiment, 
it occurs between the STM tip and the Cu-O layer when the 
tip-sample separation is small. 

In the low-bias STM image of sample A (Fig. 2(a)), we 
can identify a bond-oriented, 2-d superstructure through- 
out the entire STM image. The 2-d superstructure ap- 
peared with the same pattern in both STM measure- 
ments at positive and negative biases. Figure 2(b) shows 
the line profiles of STM images, taken along the solid 
line shown in Fig. 2(a), for various bias voltages. The 
2-d superstructure with a period of 4a appears clearly 
below V s ~ 100 mV in addition to the underlying prim- 
itive lattice, and the period of 4a is almost independent 
of bias voltage V s . The superstructure is more intense at 
lower biases, while it becomes very weak above V s = 100 
mV. Part of the low-bias STM image, taken on sample 
B at 14 = 10 mV, is shown in Fig. 3(a). The 2-d su- 
perstructure appears locally on a nanometer scale, not 
throughout the entire STM image. In Fig. 3(b), a line 
profile of the STM image is shown for the area where the 
clear 2-d superstructure appears locally and compared 
with that of sample A. This line profile shows that the 
period of the superstructure is also 4a but its amplitude 
is much smaller than that observed for sample A. 

The 2-d superstructure can also be confirmed in the 
Fourier map F(q x ,q y ) of the low-bias STM images. The 
Fourier map F(q x , q y ) of the STM image, taken on sam- 
ple A at V s = 30 mV, shows that the main Fourier 
peaks associated with the 2-d superstructure appear at 
q = (1/4,0) and (0,1/4), as shown in Fig. 4(a). This 
means that the period of the 2-d superstructure is 4a x 4a. 
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FIG. 2: (Color online) (a) Part of a low-bias STM image of 
sample A, measured at a bias voltage of V s — 30 mV and an 
initial tunneling current of It = 0.08 nA at T ~ 9 K, showing 
a 4a x 4a superstructure, together with individual atoms. The 
inset is part of a high-bias STM image of sample A, measured 
at V s = 600 mV and I t = 0.3 nA at T ~ 9 K. The image 
shows a one-dimensional (1-d) superlattice, inherent in the 
Bi-O layer, with missing atom rows, (b) Line profiles taken 
along the solid line in the STM image (Fig. 2(a)) for various 
bias voltages. The solid line is cut perpendicular to the b 
axis, that is, 45 degrees from the orientation of the 4a x 4a 
superstructure so that the 1-d superlattice of the Bi-O layer 
does not obscure the profile of the 4a x 4a superstructure. 
Note that, in the line profile for the lowest bias, the spatial 
variation due to the underlying host lattice is partly cut over 
the intense 4a x 4a superstructure because of saturation of 
the STM amplifier. 



The Fourier spot at q = (1/4, 0) is stronger than the spot 
at q = (0, 1/4), indicating that the 2-d superstructure is 
anisotropic. In Fig. 4(a), a line cut of the Fourier map 
along the (jr, 0) direction is also shown for sample A as a 
function of the bias voltage V s . Each line profile is nor- 
malized with the intensity of the Bragg peak at q = (1, 0). 
Weak Fourier peaks also appear at q = (3/4,0) in ad- 
dition to the strong main peak at q = (1/4,0). Both 
q = (1/4,0), (3/4,0) peaks are most intense at the low- 
est bias (20 mV), but they decrease rapidly with V s and 
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FIG. 3: (Color online) (a) Part of a low-bias STM image of 
sample B, measured at V s = 10 mV and It = 0.08 nA at T ~ 9 
K. The inset is part of a high-bias STM image of sample B, 
measured at V s = 600 mV and It = 0.3 nA at T ~ 9 K. Note 
that it shows almost no missing atom rows in contrast with 
that of sample A (the inset of Fig. 2(a)). (b) Line profile of 
STM image taken along the solid line in Fig. 3(a) at V s — 10 
mV (solid line). For comparison, the line profile at V B = 30 
mV for sample A (Fig. 2(b)) is also shown (dashed line). The 
dashed line for sample B is a guide to the eye. 



become very weak above V s ~ 100 mV. It should be noted 
here that these Fourier peaks show no change in position 
and no broadening as V s increases, providing evidence 
that the present 4a x 4a superstructure is nondispersive. 
The 4a x 4a superstructure can be observed even above 
T c (= 72 K) in sample A, as shown in the inset of Fig. 
4(a). 

In Fig. 4(b), the Fourier map of the STM image and its 
line cut along the (tt, 0) direction are shown for sample 
B. The q — (1/4,0) Fourier peak appears up to V s = 50 
mV, with a very weak peak at q = (3/4,0). The inten- 
sity of the q = (1/4,0) peak, normalized with the Bragg 
peak intensity at q — (1,0), is much weaker than that 
of sample A. The q = (1/4,0) peak decreases with V s 
and becomes very weak above V s ~ 50 mV. In the line 
cut of the Fourier map for sample B as well as sample 
A, peak structures are observed at q < 0.2. However, 
the intensity of these peaks is almost independent of V a , 
meaning that these structures are irrelevant to the 2-d 
superstructure we focused on (Figs. 4(a) and (b)). The 
present nondispersive 4a x 4a superstructure observed in 
samples A and B is essentially the same as the nondis- 
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FIG. 4: (Color) (a) 2-d Fourier maps of the STM images 
of sample A, taken at 14 = 30 mV at T ~ 9 K (the lower 
inset) and ~ 82 K (> T c ) (the upper inset), and cut along the 
(0, 0)-(7r, 0) line in the Fourier maps at various bias voltages. 
The Fourier amplitude is normalized by the intensity of the 
Bragg peak except the amplitude at the lowest bias, which is 
normalized so that its background level agrees with those for 
other biases, (b) Fourier map of the STM image of sample B, 
taken at V a — 30 mV at T ~ 9 K (the inset), and cut along 
the (0, 0) - (it, 0) line at various bias voltages. The Fourier 
amplitude is also normalized by the intensity of the Bragg 
peak. 



persive 4a x 4a charge order with the internal structure 
of 4a/3 x 4a/3 reported by Hanaguri et al. for lightly 
doped Na-CCOC. 1 

Shown in Fig. 5(a) is part of a low-bias STM image 
taken on sample C at V s = 30 mV, which was cut from 
the single crystal /3 (p ~ 0.13, T c ~ 78 K). The 2-d 
superstructure also appears throughout the entire STM 
image, as in sample A. Figure 5(b) shows the line profiles 
of STM images taken at bias voltages of 30 mV and 80 
mV. The superstructure with a period of ~ 4a clearly 
appears at 30 mV, but its amplitude is much weaker than 
that observed for sample A. No superstructure appears 
at a bias of 80 mV. On the other hand, it is difficult to 
identify 2-d superstructure in the low-bias STM image 
of sample D, which was cut from the same single crystal 
(/?) as sample C (Fig. 5(c)). This means that the 4a x 4a 




FIG. 5: (Color online) (a) Part of a low-bias STM image of 
sample C, measured at K = 30 mV and It = 0.09 nA at T ~ 9 
K, showing a 4a x 4a superstructure. The inset is part of a 
high- bias STM image of sample C, measured at V s = 300 mV 
and It = 0.3 nA at T ~ 9 K. (b) Line profiles, taken along the 
solid line in STM image (Fig. 5(a)), at different bias voltages 
(the solid lines). For comparison, the line profile at V s = 30 
mV for sample A (Fig. 2(b)) is also shown (the dashed line). 
The dashed line for sample C is a guide to the eye. (c) Part of 
a low-bias STM image of sample D, measured at V s = 30 mV 
and It ~ 0.08 nA at T ~ 9 K, showing a very weak 4a x 4a 
superstructure. The inset is part of a high-bias STM image 
of sample D, measured at V s = 300 mV and It = 0.3 nA at 
T~9K. 



superstructure is very weak in sample D. In Fig. 6, low- 
bias STM images are also shown for samples E and F, 
which were cut from the single crystal 7 (p ~ 0.14, T c ~ 
81 K). The 2-d superstructure appears clearly throughout 
the STM image of sample E, while it is rather weak in 
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the STM image of sample F. 

(a) sample E — 2 nm 




FIG. 6: (Color online) (a) Part of a low-bias STM image of 
sample E, measured at V s = 20 mV and It = 0.08 nA at 
T ~ 9 K, showing a 4a x 4a superstructure. The inset is part 
of a high-bias STM image of sample E, measured at V s — 800 
mV and It = 0.3 nA at T ~ 9 K. (b) Part of a low-bias STM 
image of sample F, measured at V s = 20 mV and I t = 0.08 nA 
at T ~ 9 K, showing a very weak 4a x 4a superstructure. The 
inset is part of a high-bias STM image of sample F, measured 
at V B = 800 mV and J t = 0.3 nA at T ~ 9 K. 



Shown in the insets of Figs. 7(a) and (c) are the Fourier 
maps F{q Xl q y ) of STM images taken on samples C and E 
at Vg = 30 and 20 mV respectively. In the Fourier map of 
sample C, the Fourier transform was carried out except 
top-right corner of the STM image measured over the 
area of 38 nm x 38 nm, where the main Fourier spot of 
the 2-d superstructure is split into a few spots because of 
local distortion of the superstructure. The main Fourier 
spots of the 2-d superstructure appear clearly in samples 
C and E, as expected from their STM images. The line 
cuts of the Fourier maps along the (tt, 0) direction are 
shown for samples C and E as a function of V s in Figs. 
7(a) and (c). The Fourier peaks associated with the 2- 
d superstructure appear at q = (~ 0.24 ± 0.01,0), (~ 
0.74±0.02,0), (0,~0.24±0.01), (0, 0.7±0.02) for sample 
C and q = (~ 0.26 ± 0.01, 0), (~ 0.74 ± 0.02, 0), (0,~ 
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FIG. 7: (Color) Line cuts of 2-d Fourier maps of STM images 
along the (0, 0)-(0, n) lines at various bias voltages for samples 
C, D, E and F. Here the Fourier amplitude is normalized by 
the intensity of the Bragg peak. The insets of (a), (b), (c) 
and (d) are 2-d Fourier maps of STM images shown in Figs. 
5(a), 5(c), 6(a) and 6(b) respectively. 



0.24 ± 0.01), (0, - 0.76 ± 0.02) for sample E. The period 
of the 2-d superstructure of samples C and E is 4a x 4a 
within experimental error, although we can not rule out 
the possibility that the superstructure is incommensurate 
in these samples. 

The 4a x 4a superstructures of samples C and E are 
evidently nondispersive, like those of samples A and B, 
as seen in Figs. 7(a) and (c). It should be noted here that 
the intensities of the main Fourier peaks associated with 
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Sample bias (mV) 

FIG. 8: (a) Averaged STS spectra of samples A, C, D and E 
on the positive bias side (V s > 0). (b) Energy (bias) depen- 
dence of the Fourier amplitude at q ~ 1/4 for sample pairs 
(A, B), (C, D), (E, F) and sample D'. The data referred to 
as samples D and D' were measured on different areas of the 
same cleaved surface. The arrow shows the sample bias V B 
corresponding to gap size A in the averaged STS spectrum 
shown in Fig. 8(a). 



the 4a x 4a superstructures of samples C and E show bias 
voltage dependence different from that of sample A; the 
peak-intensity is very weak at the lowest bias (V s = 10 
mV) but rapidly increases with the increase of V s and 
reaches the maximum at around V s = 20 — 30 mV (Fig. 
8), where the Fourier peak becomes much stronger than 
the Bragg peak (Fig. 7(a) and (c)). The intensity of the 
Fourier peak decreases above V s ~ 30 mV and becomes 
very weak above V s — 50 — 60 mV. In Fig. 7(b), the 
Fourier map of the STM image and its line cuts along 
the (w, 0) direction are shown for sample D. In these line 
cuts, we can identify the Fourier peak corresponding to 
the 4a x 4a superstructure only at V s — 30 mV, although 
it is very weak. In Fig. 7(d), the intensity of the main 
Fourier peak of the 4a x 4a superstructure is also shown 
for sample F. 

There exists the possibility that the 2-d superstruc- 
tures of samples C and E will be incommensurate, as 
mentioned above. The incommensurate superstructure 



reminds us of the weakly dispersive ~ 4a x 4a struc- 
ture of the LDOS maps which results from the SC quasi- 
particle scattering interference. 6,7 However, the present 
bias dependences of the wave numbers q = (<~ 0.24,0), 
(0, - 0.24) for sample C and q = (~ 0.26, 0), (0, - 0.24) 
for sample E are too small to be explained in terms of 
the SC quasiparticle scattering interference (Fig. 8(b)). 7 
Furthermore, it can hardly be understood in terms of 
the SC quasiparticle scattering interference that the main 
Fourier peaks of 4a x 4a superstructure are much stronger 
than the Bragg peak at V s = 20 — 30 mV. 



B. Results of STS; superconducting gap structure 
and 4a X 4a charge order 

Shown in Fig. 9(a) is the spatial dependence of STS 
spectra for sample A, which exhibits an intense 4a x 4a 
charge order throughout its entire low-bias STM image. 
Many of the STS spectra show the asymmetric V-shaped 
ZTPG, but some show a symmetric V-shaped gap with 
no peaks at the gap edge. Thus, the gap structure of sam- 
ple A is spatially heterogeneous and inhomogencous. In 
Fig. 9(b), the spatially averaged spectrum over a distance 
of ~ 35 nm is shown together with the ZTPGs reported 
for lightly doped Na-CCOC and Bi2212. 1 - 2 The averaged 
gap structure is very similar to the ZTPGs of Na-CCOC 
and Bi2212. Width of the averaged gap A , defined as 
the width between a shoulder on the positive bias side 
and zero bias V s — 0, is ~ 80 meV. In Fig. 10, the spatial 
dependence of STS spectra is shown for sample B, whose 
low-bias STM image exhibits a weak 4a x 4a superstruc- 
ture locally. Interestingly the STS spectra of sample B 
exhibit a homogeneous gap structure of the d-wave type. 
Gap width Ao , defined as half of the peak-to-peak width, 
is <~ 56 meV although it tends to be slightly enhanced 
over the region where the 4a x 4a superstructure appears 
clearly. 

In Fig. 11(a), the STS spectra are shown for sample C, 
which exhibits an intense 4a x 4a superstructure through- 
out its entire low-bias STM image. It should be stressed 
that the STS spectra are spatially inhomogeneous. Rep- 
resentative gap structures in STS spectra of Fig. 11(a) 
are shown in Fig. 11(b); the gap structure ranges from 
a typical d-wave type to an asymmetric V-shaped type 
with no peaks at the gap edges, and a gap with larger 
width tends to be accompanied by a subgap. The vari- 
ation of gap structure shown in Fig. 11(b) is very sim- 
ilar to that reported by McElroy et al. for underdoped 
Bi2212 crystals. 2 It should be noted that the gap struc- 
ture around the bottom is almost the same among all 
the STS spectra, although the entire gap structure dif- 
fers among them, as seen in Fig. 11(b). This means that 
the quasiparticle states around the nodes of d-wave gap, 
which dominate the gap structure around the bottom, 
are homogeneous, and so the inhomogeneity of the gap 
structure should be attributable to the nature of quasi- 
particle states away from the nodes, namely, around the 
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FIG. 9: (a) Spatial dependence of STS spectra at T ~ 9 K for 
sample A. Solid and dashed lines represent the asymmetric V- 
shaped ZTPG and the symmetric V-shaped gap with no peaks 
at the gap edge, (b) STS spectrum averaged over a distance of 
~ 35 nm on the cleaved surface of sample A at T ~ 9 K. Typ- 
ical ZTPG spectra of lightly doped Na-CCOC (dotted line) 
and Bi2212 (dashed line) are also shown for comparison. 1,2 



antinodes (Fig. 12). In sharp contrast to sample C, the 
gap structure of sample D, exhibiting a very weak 4a x 4a 
superstructure, is spatially homogeneous and of a typical 
d-wave type with A ~ 48 meV, as shown in Fig. 13(a). 
The STS data as well as STM data for sample D were 
obtained over two different areas of the same cleaved sur- 
face, but no different features appear in both STS and 
STM data (Fig. 8). In Fig. 13(b), the spatially averaged 
STS data of samples C and D are shown for comparison. 
Note that the gap structure around the bottom, \V S \ < 20 
mV, is the same between samples C and D, which were 
both cut from the same single crystal (3. The agreement 
means that the pairing gap structure around the node 




Sample bias (mV) 

FIG. 10: Spatial dependence of the STS spectra of sample 
B, taken at T ~ 9 K along the dashed line with a length of 
~ 30 nm in Fig. 3(a). Two-headed arrows beside the spectra 
indicate the regions where the 4a x 4a superstructure is clearly 
observed. 



is almost the same in the sample pair C and D, indicat- 
ing that the doping level is not so different between both 
samples. 

Figure 14 shows the STS spectra for sample E, which 
exhibits a clear 4a x 4a superstructure throughout its 
low-bias STM images. The STS spectra are also inho- 
mogeneous spatially, as in sample C. In Fig. 14(b) rep- 
resentative STS spectra are shown for sample E. On the 
other hand, the STS spectra of sample F, exhibiting a 
very weak 4a x 4a superstructure, is homogeneous as in 
sample D. 

It should be emphasized here that STS spectra of sam- 
ples A, C and E, which exhibit intense 4a x 4a superstruc- 
tures throughout their low-bias STM images, are spa- 
tially inhomogcncous (heterogeneous). As seen in Figs. 
8(b), 11(b) and 14(b), sample C exhibits more intense 
4a x 4a superstructure than sample E, and the STS spec- 
tra of sample C show more various types of the pairing 
gap than those of sample E, implying that the sample ex- 
hibiting more intense 4a x 4a superstructure shows more 
inhomogeneous gap structure. On the other hand, STS 
spectra of samples B, D and F, exhibiting weak and/or lo- 
cal superstructures, are rather homogeneous. Similar re- 
lation between the inhomogeneous gap structure and the 
development of the 4a x 4a superstructure was confirmed 
on another sample pair, different from sample pairs (C, 
D) and (E, F). These facts will lead us to the possibility 
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FIG. 11: (Color) Spatial dependence of STS spectra of sample 
C, taken along the dashed line with a length of ~ 30 nm at 
T ~ 9 K in Fig. 5(a). Representative gap structures are col- 
ored, (b) Representative gap structures shown by the colored 
lines in Fig. 11(a). The dashed line shows the STS spectrum 
averaged over all the STS spectra in Fig. 11(a). 



that inhomogeneous gap structure will be related with 
the development of the 4a x 4a superstructure. 

In the present study, the STS measurements were car- 
ried out on the same area of the cleaved surface where 
STM images were taken. Thus, the inhomogeneous gap 
structure, relating with the nature of quasiparticle states 
around the antinodes, will be intrinsic to the 4a x 4a 
superstructure state; namely, in the 4a x 4a superstruc- 
ture state the quasiparticle states around the antinodes 
will be modified inhomogeneously. McElroy et al. have 
reported that the nondispersive 4a x 4a superstructure, 
caused by a charge order, brings about a severe decoher- 
ence effect on quasiparticle states around the antinodes. 
They claimed that the observation of the charge order 
would be restricted to outside of the pairing gap with 
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FIG. 12: (Color online) Illustration for "the Fermi arc." Note 
that the pairing gap is inhomogeneous around the antinodes 
near (tt, 0) and (0, n). 
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FIG. 13: (a) Spatial dependence of STS spectra of sample D, 
taken along the dashed line with a length of ~ 40 nm at T ~ 9 
K in Fig. 5(c). (b) STS spectra averaged over all spectra for 
samples C and D. 



9 




■100 -50 50 

Sample bias (mV) 



FIG. 14: (Color) (a) Spatial dependence of STS spectra of 
sample E, taken along the dashed line with a length of ~ 30 
nm at T ~ 9 K in Fig. 6(a). Representative gap structures 
are colored, (b) Representative gap structures shown by the 
colored lines in Fig. 14(a). The dashed line shows the STS 
spectrum averaged over all the STS spectra in Fig. 14(a). 
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FIG. 15: (a) Spatial dependence of STS spectra of sample F, 
taken along the dashed line with a length of ~ 15 nm at T ~ 9 
K in Fig. 6(b). (b) STS spectra averaged over all spectra for 
samples E and F. 



Ao > 65 meV. 2 However, this is not the present case. 
Figure 8 shows that the 4a x 4a superstructure appears 
conspicuously within the pairing gap with Ao < 65 meV. 
If the 4a x 4a superstructure is due to some kind of lattice 
distortion, it is difficult to explain why the superstructure 
appears within a limited bias (energy) range, especially, 
associated with the gap size Ao; V s < Ao/e. Such a 
limitation of V s in observing the 4a x 4a superstructure 
means that the superstructure is electronic in origin, that 
is, due to an electronic charge order, as has been claimed 
in many preceding studies. 1-3 ' 9-11 The appearance of the 
4a x 4a charge order within the pairing gap implies that 
quasiparticles of the SC state and/or hole pairs will take 
part in causing the charge order. 

Here we pay attention to the spatial resolution of 
STM/STS tip-scanning in order to discuss the interre- 
lation between the development of 4a x 4a charge or- 
der and spatial inhomogeneity of the pairing gap. When 
the resolution of STM tip-scanning depends on wave- 
vector and happens to be very poor at q x , q y < 1/4, the 
4a x 4a charge order observed in STM experiment will be 
smeared and seemingly weakened, although the underly- 



ing lattice (each atom) can be observed clearly. In such 
a case, even if the pairing gap is spatially inhomogeneous 
in itself, it will appear homogeneous in STS experiments. 
Because the gap structure will be averaged over a wide 
area on account of poor spatial resolution of STS tip- 
scanning. Such a situation might be in samples B, D and 
F, whose 4a x 4a charge order is weak and gap struc- 
ture is homogeneous. Then, in order to examine whether 
the spatial resolution of tip-scanning is high enough for 
observing the 4a x 4a charge order in the present study, 
we focus on the 1-d superlattice, whose wave vector Q 
points to (1, 1) direction with \Q\ ~ l/(\/2 x 5). This is 
because 2-d charge order's wave- vectors q = (~ 1/4,0), 
(0, ~ 1/4) have a similar component (|qq| ~ X/{y2 x 4)) 
to the 1-d superlattice's wave-vector (|Q| ~ l/(\/2 x 5)) 
in the (1, Indirection. As shown in Fig. 16, the FT spots 
associated with the 1-d superlattice clearly appear not 
only in sample C exhibiting intense 4a x 4a charge or- 
ders but also in sample D exhibiting weak 4a x 4a charge 
order. Similar results were also obtained for other sam- 
ple pairs (A, B) and (E, F). These results indicate that 
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there is no essential difference in the spatial resolution 
of STM/STS tip-scanning (Q and q) among samples A 
to F. Hence the 4a x 4a superstructures in samples B, D 
and F will not be weakened by poor spatial resolution of 
the tip-scanning but intrinsically weak . 

We can also examine the spatial resolution of the tip- 
scanning from the spatial dependence of the pairing gap. 
As mentioned above, if spatial resolution of STM tip- 
scanning is very poor, the spatial dependence of the gap is 
seemingly homogeneous in STS experiments, even though 
the gap structure is intrinsically inhomogeneous. In such 
a case, the observed homogeneous gap would be essen- 
tially the same as the spatially averaged one of inhomoge- 
neous gap data, which could be obtained in STS measure- 
ments with good spatial resolution of STS tip-scanning. 
To check this point, we compare the the spatially aver- 
aged gap obtained for sample C with the averaged gap 
for sample D, because the doping levels of both samples 
are almost the same, as mentioned above. However, Fig. 
13(b) shows that the peak structure of the averaged gap 
is quite different between samples C and D; the peaks of 
sample D are evidently higher and sharper than those of 
sample C. The averaged gaps of samples E and F also 
show a similar tendency (Fig. 15(b)). These indicate 
that the homogeneous gap structures of samples D and 
F will be intrinsic, not caused by poor spatial resolution 
of STM/STS tip-scanning. 



sample D 




q< i,o)-(o. i) 



FIG. 16: Line cuts of the Fourier spots associated with the 

1- d superlattice in samples C and D, obtained from the insets 
of Figs. 7(a) and (b). Here the peak-intensities for samples C 
and D were normalized with their Bragg-peak intensities. The 
line cuts were taken along (l,0)-(0, 1) lines around the Bragg 
spots, as shown by dotted boxes in the inset. The inset is the 

2- d Fourier map for sample D shown in Fig. 7(a), although 
its contrast is enhanced so that the Fourier spots associated 
with the 1-d superlattice will clearly be visible in the inset 



C. Pinned 4a X 4a charge order and 
inhomogeneous gap structure 

In high-T c cuprates, since gap width Ao largely de- 
pends on doping level p, the value of Aq measured by 



STS provides information about doping level p of Cu- 
O layers near the cleaved surface. 23 ' 24 Gap width Ao of 
sample A, exhibiting the intense 4a x 4a charge order, is 
larger than that of sample B, exhibiting the weak 4a x 4a 
charge order locally. This means that surface doping level 
p of sample A is lower than that of sample B, although 
both samples were cut from the same single crystal a. 
Furthermore, low-bias STM imaging of slightly under- 
doped Bi2212 with A ~ 35 meV exhibits no 4a x 4a 
charge order, as was previously reported. 21 These results 
indicate that low doping tends to favor the development 
of the 4a x 4a charge order, at least, in the present study. 

Doping level p of sample B is lower than that of sample 
C, because A of sample B is larger than that of sample 
C (Figs. 10 and 11). Therefore, a more intense ~ 4a x 4a 
charge order could be expected to appear in sample B, 
compared with the charge order of sample C. However, 
the 4a x 4a charge order of sample B is weak and only ap- 
pears locally in the low-bias STM image, whereas sample 
C exhibits an intense charge order throughout its entire 
low-bias STM image (Figs. 3(a) and 5(a)). In addition 
to this result, the 4a x 4a charge order is much weaker 
in sample D than in sample C, although the doping level 
is not so different between the samples, as mentioned 
above. These results indicate the possibility that doping 
level will not be the only crucial factor necessary for the 
development of the nondispersive 4a x 4a charge order; 
there will be some other important factors in addition to 
the doping level. It should be emphasized here that sam- 
ple B shows homogeneous STS spectra with a d-wave gap, 
as mentioned above. The specific, homogeneous d-wave 
gap in sample B means that doping level p is rather homo- 
geneous and hole-pairs are uniformly formed throughout 
Cu-0 layers in this sample. These results suggest that 
the 4a x 4a electronic charge order will develop dynami- 
cally throughout Cu-0 layers, and it will be pinned down 
locally over the region with effective pinning centers in 
sample B. The pinning of the dynamical charge order will 
enable us to observe it in STM measurements. From the 
standpoint of this pinning picture, the marked difference 
of the 4a x 4a charge order between samples C and D, 
with similar doping levels, can be explained as the differ- 
ence in the density and/or strength of pinning centers. 

The Bi2212 crystals used in the present study belong 
to the pseudogap regime. In the pseudogap regime, the 
Fermi surface can be classified into coherent and inco- 
herent parts; the former is centered at the nodal point 
of the d-wave gap and often referred to as "the Fermi 
arc" , whereas the latter is around the antinodes, that is, 
outside the Fermi arcs (Fig. 12). 25-32 This heterogeneous 
structure of the Fermi surface in the pseudogap regime 
can provide a possible reason why parts of quasiparticlc 
and/or hole-pair states become inhomogeneous in the in- 
tense, pinned 4a x 4a charge order state; the incoherent 
electronic states around the antinodes, where the pseu- 
dogap develops at T > T c , are easily modified by external 
perturbation caused by the randomness associated with 
pinning potential of the charge order. 
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From the standpoint of the present pinning picture, the 
4a x 4a charge order can be expected to appear in STM 
images when the incoherent quasiparticle states around 
the antinodes (outside the Fermi arcs) contribute to the 
STM tunneling. In samples with Fermi arcs of a finite 
size, namely moderately underdoped samples, no 4a x 4a 
charge order will appear in the STM images at very low 
biases, where only coherent quasiparticles with very low 
excitation energy on the Fermi arcs contribute to the 
STM tunneling (Fig. 8). However, the 4a x 4a charge or- 
der will appear at higher biases, where incoherent quasi- 
particles with high excitation energies outside the Fermi 
arcs contribute to the STM tunneling, as observed ex- 
perimentally in samples C and E (Fig. 8). On the other 
hand, in heavily underdoped samples A and B with very 
tiny Fermi arcs, it is plausible that the 4a x 4a charge or- 
der appears in STM images even at very low biases. This 
is because there are incoherent quasiparticle states with 
very low excitation energies just outside the tiny Fermi 
arcs, and they contribute to the STM tunneling even at 
very low biases. 

IV. SUMMARY 

We performed low-bias STM imaging on underdoped 
SC Bi2212 crystals, and confirmed that the nondispersive 
4a x 4a electronic charge order appears within the pairing 
gap at T < T c . The present nondispersive charge order 
is consistent with the findings in the LDOS maps for the 
SC state of Bi2212 by Howald et al. and for the pseu- 
dogap state (T > T c ) by Vershinin et a/. 3 - 9 - 10 Howald et 
al. have claimed that the nondispersive charge order re- 
sults from the formation of the stripe order, though this 
scenario does not so straightforwardly explain why the 
observation of the 4a x 4a charge order is restricted to 
within the pairing gap. 9 . 10 . 33-35 rp ne a pp earance f the 

charge order within the pairing gap is not inconsistent 



with the models of pair density waves, electronic super- 
solids, paired-hole Wigner crystallization, or the coexis- 
tence of multi-type SC and spin density wave. 3 6 On 
the other hand, Vershinin et al. have claimed that the 
nondispersive <~ 4a x 4a charge order at T > T c is a hid- 
den order of the electron system in the pseudogap state 
(T > T c ). 3 In that case, the observation that the nondis- 
persive charge order survives even in the SC state means 
that the hidden order of the pseudogap state will remain 
essentially unchanged down to below T c . It is urgently 
desired to elucidate how the charge order in the pseudo- 
gap state (T > T c ) evolves into the nondispersive one in 
the superconducting state. 

We pointed out the possibility that the sample depen- 
dence of the nondispersive 4a x 4a charge order can be 
understood qualitatively from the standpoint of the pin- 
ning picture, which indicates that the 4a x 4a charge order 
will be dynamical in itself and pinned down over regions 
with effective pinning centers. The dynamical 4a x 4a 
charge order is a possible candidate for the hidden order 
in the pseudogap regime of pure bulk crystals with no 
effective pinning centers. We also pointed out that the 
pairing gap of samples exhibiting more intense, pinned 
4a x 4a charge order is spatially more inhomogeneous. 
The inhomogeneous gap structure can be attributable 
to incoherent electronic (quasiparticle) states around the 
antinodes, where the pseudogap develops at T > T c . The 
electronic (quasiparticle) states will be largely modified 
there by randomness associated with the pinning poten- 
tial of the 4a x 4a charge order. 
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